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I. INTRODUCTION 

Relativistic mean field (RMF) theory is a powerful ap- 
proach in describing the properties of infinite nuclear matter 
and finite nuclei simultaneously lljl OS 
OdSElSlElIIimi MM- While RMF is very 
successful at around the saturation density and in finite nu- 
clei within the range of the model parameters, there are un- 
certainties in the stiffness and the high density behavior of 
the equation of state (EOS). Stiffness of EOS can be exp eri- 
mentally probed in heavy-ion collisions IpTl I22I l23l l24l l25ll 
and the giant monopole resonances ll26l l27ll . and is crucial 
in predicting compact astrophysical phenomena. Therefore, 
it is very important to develop relativistic models with con- 
straints to this uncertainty by some symmetries, and to in- 
vestigate bulk properties of nuclear systems. One of the 
most important symmetries in hadron physics is the chi- 
ral symmetry. Chiral symmetry is a fundamental symmetry 
in QCD with massless quarks, and its spontaneous break- 
ing generates hadron masses through the chiral condensate 
(qq) Ir28ll29ll30ll . which is considered to be partially restored 
in dense matter. Thus, theories of quarks and hadrons should 
respect this essential symmetry in inquiring dense hadronic 
matter, and several attempts have been made by including 
chiral symmetry in the relativistic nuclear man y-body theo- 

ries eiliiiaiillillilHStiltnltiliaa. 

The first RMF model was proposed to deal with the prop- 
erties of nuclear matter and finite nuclei [0]. In this approach, 
the meson fields (sigma and omega) are treated as classical 
limit termed as mean field approximation. The RMF model is 
extended by introducing the isovector-vector meson p and the 
non-linear self coupling terms of mesons to obtain better de- 
scriptions of nuclear matter and finite nuclei 100, 00]. The 
significant point of this extension is that its non-linear terms 
simulate the three body forces which is essential to reproduce 
the nuclear matter saturation properties in non-relativistic cal- 
culations. There are several problems in these RMF models 
to be regarded as finite density hadronic field theories. First, 
these models do not possess chiral symmetry. While the in- 
medium nucleon mass is shifted by the a field and the scalar 
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a field may be related to the chiral condensate, the self-energy 
of a is not chirally invariant. Secondly, the vacuum is not sta- 
ble with respect to the variation of a in many of the RMF 
models. While we do not have practical problems in describ- 
ing nuclear matter and finite nuclei in the mean field treatment 
where the fluctuation of the meson fields are omitted, the vac- 
uum instability is conceptually problematic. 

The effective chiral model is analogous to the RMF model. 
We start from the cj) 4 theory, in which the spontaneous symme- 
try breaking is included, and vector mesons are introduced in 
order to describe the repulsive potential at high density. When 
we naively introduce the vector meson field into the <fi 4 the- 
ory, however, it is known that the normal vacuum jumps to 
a chiral restored abnormal vacuum (Lee-Wick vacuum) be- 
low the saturation density |0[3j]], and this problem is referred 
to as the chiral collapse problem ll32tl . One of the prescrip- 
tions to avoid this problem is to incorporate a logarithmic 
term of a |0 [H [II d E [M 13 in the chiral po- 
tential (energy density as a function of a at zero baryon den- 
sity). The logarithmic a potential is first introduced to sim- 
ulate the scale anomaly in QCD, which is represented by the 
glueball dilaton field (y) which couples with the logarithm 
of a as oc - X 4 log ° 0, M M d. It is also pos- 
sible to derive it from the st rong coupling limit (SCL) of 
lattice QCD HI [3 [H [H Hj. The logarithmic poten- 
tial term, oc — log er, generally prevents the normal vacuum 
from collapsing and hence it has no instabilities. Also this 
model can describe well the even-even finite nuclei and infi- 
nite nuclear matter properties by including the vector mesons, 
their linear c oupling s to nucleons, and a self-interaction term, 

(lu^) 2 mElQlEi. InRef - C2l' for example, Schramm 
applied a chiral SU(3) RMF model with a logarithmic a po- 
tential and the glueball interaction to nuclei over the whole 
known range assuming axial symmetry, and it was demon- 
strated that binding energies of spherical and deformed nu- 
clei are well explained with the precision of 0.1 — 1 % ac- 
curacy. While these models have met phenomenological suc- 
cesses and somewhat based on QCD, one may doubt the va- 
lidity of the logarithmic potential, since it is divergent when 
the chiral symmetry is restored, a — > 0. 

Another way to avoid the chiral collapse is to introduce 
a dynamical generation of the isoscalar-vector meson mass 
through the coupling between scalar and vector mesons |0 
[IH]. Since the vector meson becomes light when the chiral 
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symmetry is partially restored, repulsive effects from the vec- 
tor meson become strong and we can avoid the chiral collapse. 
One of the problems in this theoretical treatment is the unre- 
alistically high incompressibility value, K > 700 MeV. In 
order to make moderate value of incompressibility, Sahu et 
al. introduced the higher order terms of scalar meson, a 6 and 
<7 8 lfl7l n~8l [l9ll . In this way, we can reproduce the empirical 
values of the saturation density, binding energy, and incom- 
pressibility in symmetric nuclear matter. The advantage of 
higher order terms in the chiral potential is that we have a free- 
dom to adopt weaker repulsive vector force in the nuclear in- 
teraction and therefore, one can have a choice of desirable val- 
ues of incompressibility. In the earlier works lfl7l[T8l[l9il20ll . 
the above model with higher order terms in scalar-field inter- 
actions was extensively used in dense matter as well as hot nu- 
clear matter. In all these works, the vacuum stability at large a 
values was not critically examined for all sets of parameters. 
It has been recently pointed out |Tl3i Hi that one of the pa- 
rameters in the previous work iflal has instability at large a. 
Though several sets of parameters of the above model were 
tabulated, we find that few of them overcome the instability 
at large a. Therefore, it motivates us to revisit this model and 
put the stringent constrains on parameters for stability with 
respect to sigma field. 



In this paper, we systematically investigate the vacuum 
stability and nuclear properties in the effective SU(2) chi- 
ral model having higher order terms in the chiral potential 
V a . The condition of the vacuum stability is elucidated in 
the parameter plane in the model. We calculate the EOS for 
symmetric nuclear matter with moderate choice of the incom- 
pressibility between 300 and 400 MeV, the effective masses 
around 0.85 of nucleon mass, under the constraint of vacuum 
stability. The parameters are chosen accordingly by constrain- 
ing the above conditions at nuclear saturation points. We also 
apply our model to finite spherical nuclei, and perform a naive 
dimensional analysis (NDA) to examine the naturalness of the 
effective Lagrangian. 



II. THE FORMALISM OF SU(2) EFFECTIVE CHIRAL 
MODEL 

The effective chiral Lagrangian which includes a dynam- 
ically generated mass of the isoscalar-vector field, oj^, that 
couples to the conserved baryonic current = can be 

written as dd, 

C =V> [i§ - g a (a + i-f 5 r ■ tt) - g^ub - g p p ■ r] i\> 

+ i(a M 7r • d»ir + d^a) - V„(a, tt) 

-- A F^F. V + \g 2 ^x 2 u^ 



(1) 



We introduce a chiral symmetric type interaction up to eighth 
order of the meson field which reads, 
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where x 2 = a 2 + tt 2 , / t is the pion decay constant, = 
d^w v - dvU^ and G M „ = d il p v - d v p p + gppP^ x p v are 
the field tensors of isoscalar- and isovector-vector mesons (ui 
and p-mesons), and ip, tt, a denote the nucleon isospin dou- 
blet, isovector-pseudoscalarpion, and the scalar fields, respec- 
tively. Coupling constants of nucleon with scalar and vector 
fields are introduced as g„, g^ and g p , respectively. We work 
in natural units where h — c = k B = 1. 

The interaction terms of the nucleon and vector meson with 
the scalar and pseudoscalar mesons generate the masses of the 
nucleon and vector meson through the spontaneous breaking 
of the chiral symmetry. The masses of the nucleon and vector 
meson in vacuum are given by, 



M N = g a f n 



(3) 



All through the paper, we work in the chiral limit (m n = 0) 
for simplicity. As long as we do not explicitly include tt me- 
son effects in the mean field approximation, the results with 
finite m T in nuclear matter and finite nuclei are found to be 
very similar to those in the chiral limit. We ignore these small 
differences, since the main aim of this paper is to elucidate the 
vacuum stability condition in effective chiral models. 

The paper is organized as follows: In Sec. [II] we present a 
brief formalism of the effective chiral model. We investigate 
the properties of nuclear matter and finite nuclei in Sec. [Ill] 
Basically, we determine the suitable parameters of the model 
to explain the saturation properties of nuclear matter as well 
the vacuum stability. Then we determine the equation of state 
and the properties of finite nuclei. We also examine the nat- 
uralness of the effective Lagrangian using naive dimensional 
analysis in Sec. [Ill] We summarize our results in Sec.lIVI 



where the vacuum expectation value of the a field is replaced 
with /tt. The coefficient C4 is related to the vacuum mass of 
<7 as 



2/ 2 



(4) 



The constant parameters C& and Cs are included in the higher- 
order self-interaction of the scalar field to describe the desir- 
able values of nuclear matter properties at saturation point. In 
this work, we consider the chiral limit, where the pion mass 
TOjr is zero. In the mean-field treatment we ignore the explicit 
role of 7r mesons. 

By adopting mean-field approximation, the equation of 
motion of fields are obtained from the chiral Lagrangian. 
This approach has been used extensively to evaluate the 
EOS JH HI HI Him in many of the theoretical models 
for high density matter. Using the mean-field ansatz in uni- 
form matter, the equations of motion for the vector fields {ui 
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and p-mesons) are solved as, 
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R = p 3 



9p 



(Pp ~ Pn) , (5) 



where Y = cr/f w is the reduction ratio of the chiral con- 
densate from its vacuum value. Proton and neutron densities 
(p p and p n ) are given as, p a = ~/[k^] 3 /6ir 2 , where 
is the Fermi momentum of the proton (a = p) or the neu- 
tron (a = n), and 7 is the spin degeneracy factor, 7 = 2. 
The baryon density is the sum of proton and neutron density, 

Pb = Pp+ Pn- 

The EOS is calculated from the diagonal components of 
the conserved total energy-momentum tensor corresponding 
to the Lagrangian together with the mean-field equation of 
motion for the fermion field and a mean-field approximation 
for the meson fields. The total energy density (e) and pressure 
(P) of the uniform many-nucleon system are given by, 



e=e N (M*) + V a + 



P=P N (M*)-V<r + 
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where M* = YM N = g a a is the effective mass of the nu- 
cleon. The free (kinetic) nucleon energy density and pressure, 
£jv and Pjy, depend on the effective mass M*, as well as on 
the nuclear density p B . 

The equilibrium value of the scalar field (c) is obtained 
from the equation of motion, and it is equivalent to the mini- 
mum energy density condition, de/da = 0. The equation of 
motion in terms of Y is given as, 



C 4 (l - Y 2 ) - C e (l - Y 2 ) 2 + C 8 (l - Y 2 ) 3 
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where p s denotes the scalar density, defined as 

nfc< ° ) M*d 3 k 
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In the previous works lll7Lll8lll9ll20ll . the u>N coupling and 
aui coupling was assumed to be the same, = g au> , and the 
pion decay constant / w was not introduced explicitly. The en- 
ergy density is represented by c CT (= g^/m 2 ),^ = gt/m^, 
and B and C, which are related to the a 6 and a 8 coefficients, 
respectively. In these works, implicitly given value is not 
necessarily the same as the observed one, / ff ~ 93 MeV. 
However, it is possible to map those parameters into the co- 
efficients in the present work by comparing, for example, the 



energy density, if we do not require the condition, g u = g auJ . 
The relation with their parameters with the present parameters 
are given by, 
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M 2 
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and g M = y/cZra^. In the later discussion, we also examine 
the stability of their Lagrangians. 



III. RESULTS AND DISCUSSIONS 

In this section we will determine the parameter sets at nu- 
clear matter saturation density. We will select the parameters 
by examining the stability with respect to sigma field. We will 
then use stable parameter sets to find the properties of finite 
spherical nuclei and then examine the naturalness of the ef- 
fective Lagrangian by performing the naive dimensional anal- 
ysis. Here, we use constants M N = 938 MeV, /„■ = 93 MeV, 
m u = 783 MeV, m p = 770 MeV and g a = M H /f„. 



A. Fixing parameters in symmetric nuclear matter 

In the effective chiral model, the chiral symmetry relates 
the interaction parameters and hadron masses, and reduces the 
number of parameters, as shown in Eq. ([3]). In the present 
treatment, we have five parameters, g u , g p ,d,Ce and C$- 
Here we determine three parameters, the nucleon coupling to 
the vector field, g u , and the coefficients in the scalar poten- 
tial terms, C4 and Cq, in symmetric nuclear matter. These 
parameters are obtained as functions of Cs and the nucleon 
effective (Landau) mass M*(po), by fitting the empirical sat- 
uration point, (po, Eo/A), where Eq/A is the binding energy 
per nucleon at saturation density, p B = po. The saturation 
point plays a decisive role on finite nuclear binding energies 
and radii, thus we will adjust them to reproduce finite nu- 
clear property. Moreover, the incompressibility K and the 
nucleon effective mass M*(po) are the keys in EOS around 
the saturation point, as well as at high densities. The nuclear 
incompressibility is somewhat uncertain at saturation point. 
The desirable values of effective mass and nuclear matter in- 
compressibility are chosen in accordance with recent heavy- 
ion collision data lf2ll I22I l23i l24l l25tl . In our calculation, we 
have examined several parameter sets corresponding to each 
incompressibility and effective mass in the range of 200 — 400 
MeV and (0.8 - 0.9)Mat jH], respectively, to observe the 
sensitivity of EOS in the high density region. 

First, we introduce and examine the vacuum stability of chi- 
ral potential as the constraint on the parameters. Vacuum sta- 
bility condition of the present effective chiral models can be 
examined as follows. The vacuum energy density is given in 
V a in Eq. ©, which is rewritten as 



V* X 
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-4 = — f(X) , f(X) = —X 2 + —X + — , (13) 
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Y(p )=M N (p )/M N 



FIG. 1: (Color online) as a function of Y(po) = M^(po)/M N 
(solid line), and Ce (dashed line), Cs (dotted line) values on the vac- 
uum stability boundary. 
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FIG. 2: (Color online) Vacuum stability region (shaded area) and the 
incompressibility K in the (Y(po) = MJJpo)/M N , Cs) plane. We 
have displayed several other models Il9ll20ll (open symbols) with 
<f> 4 theory (filled square) and present results (filled circles). 



where X = (Y 2 — 1). In stable cases, V a must be always 
positive in the range X > — 1 except for the vacuum X = 
(i.e. (7 = at which V a = 0. Provided that C 8 > 0, 
the stability is ensured when one of the following conditions 
is satisfied; (1) the discriminant of f(X), D = (C 6 /3) 2 - 
4(C 8 /4)(C 4 /2), is negative, (2) the discriminant of f(X) is 
positive, but the solution of f(X) = are in the range X < 



-1, or (3) in the case of C 8 = 0, C* 6 > and /(-l) > 0. 

By fitting the saturation point, we can fix C4, Cq and as 
functions of C 8 > and the effective mass at normal density, 
Y(p Q ) = M*/M N ~ 0.8 - 0.9. First, we give Y(p ), then 
g^ is uniquely determined. From Eqs. © and ©, we find that 
the enthalpy density is free from V„, 



P = 



m 2 R 2 



rri^Y 2 

Pb(M n -B/A) {p B =po,P P 



e N (M*) + P N (M*) (14) 
Pn) ■ (15) 



In the second line, we have used the fact that P = at the sat- 
uration density. This equation only depends on one parameter 
g w , and we can fix it from the saturation property. In Fig. Q] 
we show g^ value as a function of Y(po) = M*(po)/M N . 
The energy gain from the a meson is small for larger val- 
ues of Y(po), then the repulsive potential is also chosen to be 
small to reproduce the binding energy B/A at pq. Thus g^ is 
a decreasing function of Y(po). Next, we give the value of 
Cs > 0. For a given set of (Y(po), Cs), we can solve the con- 
dition, e(po) = Po{Mn — B/A) and de/da = at p B = p 
(Eq. (fTUll), with respect to C4 and Cq. 

From these coefficients, we examine the vacuum stability 
condition. In Fig.Q] we show Ce and Cs values on the vacuum 
stability boundary as a function of Y(po). These values are 
equivalent to the minimum Ce and Cs values for each Y(po). 
At Y(p ) = 0.781 = >04, 4> A theory is realized, i.e. C 8 = 
Ce = 0. For larger value of Y at p B = po, repulsion from w is 
chosen to be smaller as shown in Fig.Q] and larger repulsion 
in V a is required. This repulsion can be generated by negative 
Ce or positive C 8 , since (Y 2 — l) 3 and (Y 2 — l) 4 are negative 
and positive for Y < 1, respectively. Thus for small values of 
Cs, negative Ce values are required at Y(po) > Y^i, and to 
keep the vacuum stability, there exists the minimum value of 
Cs- 

In Fig. 12 we show the vacuum stability region (shaded area) 
in the (Y(po), Cs) plane. We have examined the stability of 
the parameter sets proposed in previous works ifTii [l9l I20I1 . 
We show these parameter sets by open symbols in Fig.|2]and in 
TablelU Unfortunately, most of the parameter sets which gives 
medium incompressibility (200 McV < K < 400 MeV) are 
unstable in vacuum against the variation of a, marked with 
"U" in the last column of Table U Only one parameter set 
(SIPP-V) fulfills the vacuum stability condition and gives a 
medium K value. 

Parameter sets which we propose and examine in this paper, 
STO-z (? = 1, 2, • • • 5), are tabulated in Table U and shown in 
filled circles in Fig. [2] We have chosen two values of M* /M N 
(0.835 and 0.85), and three values of C 8 (20, 40 and 60). The 
combination (M*/M N ,C 8 ) = (0.85, 20) is close to the sta- 
bility boundary, and we do not adopt them. All of these pa- 
rameters give stable chiral potentials, and the incompressibil- 
ity are in the range, 200 MeV < K < 400 MeV. We adopt 
the saturation point (p ,E /A) = (0.14 fm~ 3 ,-(14.5 — 
14.6) MeV), which is found to explain the binding energies 
of heavy nuclei reasonably well as discussed in the next sub- 
section. 

In Fig. [3j we show the chiral potential V a , in STO-5 as an 
example. We also show the chiral potential in the <fi 4 theory, 
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TABLE I: Effective chiral model parameter sets. Stability of the model is also shown; the vacuum is stable against the variation of a in models 
with "S", and the vacuum is unstable at lame a in models with "U". 
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404.300 


173.400 


706.2 


4.53 


s 



TM1 |H, SCL d H and SO 03, for comparison. The 
chiral potential in STO-5 behaves similarly to that in SO in 
the region a < f n . At larger a values, unstable parameters 
give negative chiral potentials in the region a > f n as shown 
in the SO case. In Fig. |4] we show the EOS in STO-5 in 
comparison with other EOSs. We find that EOS in STO-5 is 
reasonably soft, and comparable to those in TM1 and SCL, 
which explains the bulk properties of finite nuclei. 



B. Finite nuclei 

In describing finite nuclei, it is numerically preferable to 
represent the Lagrangian in the shifted field <f = / w — a 
and to separate the a mass term from the chiral potential V a , 
since the boundary condition is given as ip — > at r — > 00. 
In addition, it is necessary to include the photon field which 
represents the Coulomb potential. Here we take the static 
and mean-field approximation for boson fields, then RMF La- 



grangian can be written as follows, 

£™ F [id - M* N {ip) ~ i°U v (u, R, A)] i> 

- \ (v<p) 2 - \ m l<p 2 - v A<p) 

+ \{VR? + \™ 2 P R 2 + \{VA) 2 , (16) 

where 

MZ(<p) = M N - 9<r <p , M»= S L(/^) 2 , (17) 

U v (w t R,A) =g u u> + g p T 3 R+el-tp-A, (18) 

V v = V a - im> 2 . (19) 
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FIG. 3: (Color online) The chiral potential as a function of a in 
the cj> 4 theory (dotted line), SCL (dashed line) [H, SO (open 
square) 0], and STO-5 (solid line). Results of a non-chiral model, 
TM1 (filled circles) d are also shown for comparison. 



sity given as, 

- - J {-g a <pp s + gu^pB + gpRpr + e 2 Ap p B } dr 

V v -\^-?fv{U-v)u 2 }dv (25) 

where we use the Eq.(|20l>-(f24b to calculate second order 
derivatives of meson fields. Nucleon single particle states are 
specified by the radial quantum number i, isospin a(= p, n), 
and the angular momentum quantum number, n = I {k = 
-(I + 1)) for j = I - 1/2 (j = I + 1/2). The number of 
occupied nucleon is represented by n°™, which is equal to 
2\k\ = 2j + 1 for filled single particle states. We solve the 
self-consistent coupled equations (1201) -(124b by iteration until 
the convergence of total energy is achieved. In this work, we 
assume that the nuclei under consideration are spherical, then 
the nucleon wave functions are expanded in spherical har- 
monic basis as follows, 



> 

01 




P B (fm" 3 ) 



FIG. 4: (Color online) Energy per nucleon as a function of the baryon 
density. Meaning of the lines and symbols are the same as in Fig. [3] 



The field equations of motion derived from this Lagrangian 
read, 



-ia ■ v + 0M* + U v ] i\) = eiip , 



(20) 



(- A + ml) ip = g aPs - -^-gL (U-f) , (21) 

(- A + ml) u = guPs+glujip (2U-<f) u , ( 22 ) 

(-A + m 2 p )R = g p p T , (23) 

-AA = epl , (24) 





E 

iK 


/ occ 

^47rr 2 


E 

IK 


/ occ 



-WJr]<S> 



fx I I I 

k i n 



Co 



f )(\G? K (r)\ 2 + \F? K (r)\ 2 ) , 



? = E US )(\G? K (r)\ 2 -\F? K (r)\ 2 ) , 



(26) 
(27) 
(28) 



where ( a represents the isospin wave function, a = p, n. 

In comparing the calculated results in mean-field models 
with the experimental binding energies and charge radii, we 
have to take account of several corrections. In this work, we 
consider the center-of-mass (CM) kinetic energy correction 
on the total energy, and CM and nucleon size correction on 
nuclear charge rms radius in the same way as that adopted in 
Ref. JU . The CM kinetic energy is assumed to be 



Ezpe = 



2AM 



N 



-hio^hlA- 1 / 3 , 
4 4 



(29) 



where Pcm = J2i Pi i s tr, e CM momentum. This correc- 
tion gives an exact result when the state is represented by a 
harmonic-oscillator wave function, and we assume that it also 
applies to the RMF wave functions. The CM correction on the 
proton rms radius is written as 



6(r 2 p ) = -2(R CM • 
3h 



2AM n uj 
2<^> , (ft 



A 



A 



(for heavy nuclei) , 
(for light nuclei) , 



(30) 



where p s = p% + p% p B = p p B + p%, p T = p p B - p n B de- 
note scalar, baryon and isospin densities of nucleons, respec- 
tively. Total energy is given by the integral of the energy den- 



where R p = J2ie P r i/Z is the proton CM position, and (r 2 ) 
and (r 2 ^) represent the proton and matter mean square radii, 
respectively. We assume again that harmonic-oscillator results 
applies for heavy nuclei. For light nuclei, we evaluate the 
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FIG. 5: Neutron single particle levels of ' 
in the chiral limit, and STO models. 
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correction in RMF wave functions, and we consider only the 
direct-term contributions. The charge rms radius is obtained 
by including the finite size effects of protons and neutrons, 



>ch/ 



size/ P 



N 
~Z 



r 2 ) 

size/ 



(31) 



denotes the size of proton or neutron and is 
equal to (0.862fm) 2 and (0.336fm) 2 , respectively. We evalu- 
ate the binding energies and charge rms radii with these cor- 
rections, and the pairing energy for open-shell nuclei are ne- 
glected. 

In describing finite nuclei, the isospin dependent interaction 
is an important ingredient. In uniform nuclear matter, we can 
obtain the symmetric energy coefficient a sym by expanding 
the energy density around the symmetric nuclear matter, 



n 2 k 3 
9p K F 



*sym 



37T 2 7 



MS 



(32) 



where k F = (6ir 2 p B /j) 1/3 (p B = p p + p„,7 = 4) is the 
Fermi momentum in symmetric nuclear matter. In order to 
reproduce the empirical value of the symmetry energy coef- 
ficient, a sym = 32 ± 6 MeV l35ll . the above relation gives 
g p = 4.625 for and M* = 0.85M. In the present work, we 
have fixed g p by fitting the binding energies of heavy nuclei. 
We find that g p ~ 3.5 is appropriate. 

We have tuned the model parameters so as to explain the 
binding energies of heavy nuclei. In TableHIl we tabulate finite 
nuclei results. We also summarize the STO parameter sets and 
the resulting incompressibility in Table Q] The binding energy 
per nucleon in heavy nuclei and the incompressibility param- 
eters would be in the acceptable range. On the contrary, we 
cannot reproduce the binding energies of light nuclei simulta- 
neously with heavy nuclei. Calculated values of charge rms 
radii underestimate the data as long as we vary the po value in 
the acceptable range, po = 0.14 — 0.16 fm" 3 . 

While we show the results obtained in the chiral limit 
(TOjt = 0), the results are not modified much with finite pion 



TABLE II: Calculated results for saturation property of symmetric 
nuclear matter, B/A, charge rms radii of stable nuclei. 

Saturation property 





STO-1 


STO-2 


STO-3 


STO-4 STO-5 


Exp. 


Po 


0.14 


0.14 


0.14 


0.14 


0.14 




E 


-14.6 


-14.5 


-14.6 


-14.6 


-14.5 




K 


318.5 


376.0 


389.5 


402.3 


327.2 




B/A 


12 C 


9.38 


9.40 


9.37 


9.26 


9.30 


7.68 


16 


10.9 


10.9 


10.9 


10.7 


10.8 


7.98 


28 Si 


9.53 


9.57 


9.87 


9.47 


9.47 


8.45 


40 Ca 


9.91 


9.85 


9.56 


9.80 


9.86 


8.55 


48 Ca 


9.73 


9.74 


9.74 


9.68 


9.70 


8.67 


90 Zr 


9.13 


9.11 


9.15 


9.09 


9.11 


8.71 


lie Sn 


8.83 


8.81 


8.85 


8.80 


8.81 


8.52 


196p b 


7.87 


7.85 


7.91 


7.86 


7.86 


7.87 


208p b 


7.87 


7.85 


7.91 


7.86 


7.87 


7.87 


Charge rms radii 


12 C 


2.27 


2.27 


2.27 


2.27 


2.28 


2.46 


16q 


2.44 


2.44 


2.44 


2.44 


2.44 


2.74 


28 Si 


2.95 


2.94 


2.95 


2.95 


2.95 


3.09 


40 Ca 


3.30 


3.30 


3.30 


3.30 


3.30 


3.45 


48 Ca 


3.39 


3.39 


3.39 


3.39 


3.39 


3.45 


90 Zr 


4.21 


4.20 


4.20 


4.20 


4.21 


4.26 


116 Sn 


4.59 


4.59 


4.59 


4.59 


4.60 


4.63 


196p b 


5.48 


5.47 


5.47 


5.47 


5.48 




208p b 


5.56 


5.55 


5.54 


5.54 


5.56 


5.50 



mass in the mean field approximation. This is examined in 
an RMF model with the chiral SU(2) logarithmic potential 
(SCL2) with finite lfl3ll . In Fig. [5] we present calculated 
neutron single particle levels of 208 Pb in the chiral limit. Not 
only in the level structure but also in their energies, there is lit- 
tle differences between the SCL2 results with finite and zero 
pion masses. Therefore in the mean field approximation, we 
can safely discuss the properties of the RMF Lagrangian in 
the chiral limit. 

We also find in Fig.[5]that the present model (STO-1) show 
smaller Is splitting compared with other RMF models due to 
large effective nucleon masses. One of the original motiva- 
tions to use RMF models for nuclei was the large Is splitting 
naturally generated from the large scalar and vector poten- 
tials additively. In the recently developed chiral RMF mod- 
els U^O, Is splittings are evaluated to be smaller than em- 
pirical values. This may be suggesting the need to include ex- 
plicit pion effects. There are some discussions regarding the 
contribution to the Is splitting 136] where one-pion exchange 
tensor force and two-pion exchange with the excitation of vir- 
tual A(1232) isobars are taken into account for explaining 
Is splitting property of nucleon and hyperon simultaneously. 
The is-like roles of the tensor force or pions are discussed also 
in light nuclei l37ll . We have not taken account of these pion 
effects in our present calculations, then 7r treatment may be 
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required in order to resolve the Is splitting problem in chiral 
RMF models and effective chiral models. 



find that the meson-nucleon and old couplings are natural, but 
the self-interaction coefficients in a are not natural. 

IV. SUMMARY AND CONCLUSIONS 



C. Naive dimensional analysis 

The present STO model is a kind of effective field theory, 
contains higher order terms and is non-renormalizable. Then 
it would be valuable to examine the naturalness in the naive di- 
mensional analysis (NDA) (Mill 5(1 El 51. It is found that 
the loop contributions with the momentum cutoff A ~ 1 GeV 
generate the followin g ter ms with dimensionless coefficients 
Qmnp of order unity, [|H IM 53 51 



Ant ~ 

_£_ 

J* 



Clmnp ( 1pTlp \ 



m\n\p\ 



(33) 



where T denotes the 7 and r/2 when necessary. 

An effective theory having terms in Eq. d33l l is considered 
to hold naturalness, when all the dimensionless coefficients 
Cimnp are of order unity. In the present effective Lagrangian, 
we obtain the following dimensionless coefficients, 



— ^ = — u - 94 1 



C1010 = ^ ~ 0.56 
C 100 i = ^~0.64 



A 



0120 



^Ul 2m 2 



A 2 



A 2 



1.2, 



_ IglJl _ 2m 2 

L-0220 — — — m ~ i-l 



A 2 



A 2 



31 /4 
C0300 = -j^ 3! I — Ce — C4 I ~ — 4 

Omoo = ^4! (2C 8 - 2C e + ic 4 



0500 



f 1 

■I TV 

A 2 

f'2 



5! (-4C 8 + C 6 ) - -280 , 



C 06 oo = |§ 6! f 3C 8 - ^C 6 j ~ 1200 

C0700 = -|§ 7!C 8 ~ -2600 , 
/2 8! 

Cosoo = T9 TT^S ~ 2600 . 
A^ 8 



(34) 



We show the results in STO-5, and adopt A = 1 GcV. We 



In this paper, we have investigated the properties of nuclear 
matter and finite nuclei in the effective chiral model with a 6 
and cr 8 terms. The nucleon-vector meson coupling is found 
to be uniquely determined as a function of the effective mass 
at normal nuclear matter density, Y(po) = M*(po)/M N , and 
we have specified the region of stability in the (Y(po), Cs) 
plane, where C$ is the coefficient of the cr 8 term. We can 
find the parameter sets which satisfies the vacuum stability 
condition and results in moderate incompressibility, K = 
(200 — 400) MeV. The incompressibility is found to be domi- 
nated by the nucleon effective mass, and M*(p )/M N > 0.83 
is necessary in order to obtain moderate K, as far as the vac- 
uum stability is required. 

The obtained effective chiral model with higher order terms 
in a is applied to finite nuclei for the first time. We can explain 
the binding energies of heavy nuclei (Sn and Pb) reasonably 
well, while we overestimate the binding energies of light nu- 
clei. This may be because the nucleon-vector meson coupling 
is small, g^N ~ 6 compared with other RMF models which 
explains nuclear binding energies in a wide mass range, such 
as NL1 (g wN = 13.285) H II, NL3 = 12.868) 

TM1 ( gujN = 12.6139) H, and SCL ( gujN = 13.02) SB- 
Smaller g^N value is compensated at higher densities where 
the chiral symmetry is partially restored and to mass de- 
creases, but light nuclei are more sensitive to the EOS at lower 
densities. 

We have also performed the naive dimensional analysis 
(NDA) H HI m 51 51 of the present model. Moder- 
ate K value of around 300 MeV requires the cr 8 coefficient 
Cg > 20 as found in Fig. [2] This value corresponds to 
Cosoo ^ 870, and the model cannot hold naturalness. In order 
to construct effective chiral models having moderate incom- 
pressibility, vacuum stability and naturalness simultaneously, 
it would be necessary to introduce other types of interaction 
terms other than polynomial forms of a. Works in this di- 
rection would be valuable for the understanding of the chiral 
properties of the QCD vacuum. 
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